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Coincidence-structures in the transverse plane of Type-II spontaneous parametric down- 
conversion carrying orbital angular momentum are obtained. Azimuthal symmetry breaking around 
the pump beam direction reveals itself on these quantum images. Analytical expressions for the am- 
plitude probability of the down conversion process are shown including the nonlinear polarizability 
, components. 
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C\l ; INTRODUCTION 

^ • Entanglement of photon states carrying orbital angular momentum (0AM) is a new tool in quantum optics. These 

' ' ' states have been experimentally produced by Spontaneous Parametric Down-Conversion (SPDC) pj]. They were first 

^^O predicted 0] subjected to certain symmetry conditions that are still subject to some controversy. Experimentally, 

' these states can be created under the condition of a-posteriori or a-priori 0AM imprinting by appropriated masks or 

' filter converters. For the a-posteriori imprinting, the pump beam does not need to carry 0AM but an appropriate 
0AM mode converter I is inserted onto one of the down converted beams. The conjugate beam is expected to acquire 

4!) I an 0AM of —I. For the a-priori imprinting, the pump beam mode is set in a 0AM state, say, I. However, as discussed 

^1 in the initial 0AM I may or may not be transferred to the SPDC. Discussions about if this transfer of 0AM occurs 

^ . or not and if special conditions are needed to have 0AM transfer has populated the recent literature sometimes with 

qh' conflicting answers: "yes", "no", "under certain conditions" (See some references in Q). Sometimes, conclusions 

'— 'i derived from particular cases say, under conditions adequate for low wave vectors or paraxial cases (e.g., Q), can be 

^ ■ mistakenly taken as being general. Part of this non-uniform understanding on this subject derives from oversimplified 

^ Hamiltonian or wave-states considered (simplified "backbone" models). The difficulties to derive general conclusions 

' taking into account realistic phase matching and light-matter coupling through the non-linear polarizability tensor 

0^ . sometimes makes these time consuming tasks not appealing. Therefore, it is common that a single proportionality 

2? ' constant replaces an involved angular dependence resulting in a non realistic description. Instead of being just 

. technical details, these more complex dependences may be crucial to a full understanding of these processes including 

■ entanglement between conjugate photon pairs. Phase matching is also frequently assumed under simplified conditions 



o 
o 



not sufficient to treat in detail pump modes with amplitudes more complex than a simple gaussian intensity profile. 

Although not claiming a complete analysis of this problem, this paper does not take for granted many of these 
commonly oversimplified assumptions. It tries to discuss the problem of transfer of 0AM between the pump beam 
. ^ ' Eind the SPDC photons with reasonable detail. Phase matching conditions dependent on x'^^^ ^-re not oversimplified; 
. also detailed is the dependence of SPDC on x^"^^ ■ This should provide the reader with enough material to help 
' his own work or, at least, to show fundamental elements necessary in this trade. It may also stimulate others to 
improve our understanding on this area. Certainly, many future applications will demand a careful understanding the 
entanglement carried by these non-linear processes beyond current treatments. Even some small immediate rewards 
can be obtained such as a method to obtain the non-linear coefficients in x*^^^ by comparison of experiment and theory. 

This paper was written trying to give a newcomer to this field a straightforward view of the involved elements 
starting from the wave state amplitude describing SPDC. It describes all components involved in this amplitude and 
detail aspects of phase matching and simplifications adopted including the description of phase matching in two steps 
(longitudinal and transverse conditions). The non- linear tensor x'^^ and the resulting non- linear polarizability of the 
medium are discussed, including their transformations between crystal reference system and the laboratory reference 
system. This way, the reader can dedicated himself to the main aspects presented instead of having to work out 
technical details. Certainly, for some, most of the presented aspects may be trivial and can be skipped without 
problems. However, any disagreement between views may result in an overall enlightenment for all. 
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WAVE STATE 



The Hamiltonian describing SPDC can be easily found in the hterature (e.g., ^5i|). It describes free propagating 
photons from a pump and from the down-conversion process that occurs due to the non-hnear hght-matter interaction 
occurring in a ideally transparent medium (e.g, crystalline medium within the crystal band gap). A pump photon 
will excite the non-linear medium through a very fast interaction with virtual electrons and decay either into a similar 
pump photon or into two conjugate photons-historically called signal and idler. These virtual interactions as well 
as the propagation of the SPDC photons occur in the non-isotropic medium with specific symmetries defined by the 
crystal class involved. Therefore, medium symmetries are built explicitly in x'^^ and implicitly in x^^\ defining the 
non-linear interaction and the light propagation in the medium. The Hamiltonian for these processes usually neglect 
coupling to the lattice possibly intermediated by electrons. Although the fast interaction times indicate that this 
coupling to the lattice should be negligible, this possibility should not be ruled out in general. 

At this point, one could remind the reader the uncertainty relationship connecting angular momentum and phase 
uncertainties 

Ai,A0 > ^ (1 - 27rP(0o)) . (1) 

It states that a large uncertainty in angle allows Lz to have a small uncertainty. Consequently, in order to guarantee 
precision in Lz, uncertainty in has to be maximum. Could the medium symmetries (built in x'^' and x^^^) cause 
restrictions on the photon interactions or propagation that could somewhat constrain the associated azimuthal angle 
(j> and therefore do not allow its maximum uncertainty to be achieved? This problem has been discussed in Ref. [3| 
but it will detailed here for clarity. 

The interaction Hamiltonian Hj gives the wave state |-0 {t)) = exp {—i/fi) J^_^ Hi{T')dT' |0), from which successful 
spontaneous photon conversion in first order is 

l^(i)> A'^'A'fc^^.,.'(k,k')a^(k,s)at(k',s')|0) . (2) 



The probability amplitude for signal and idlers at (k, k') is Fg^s' (k, k') = ^k,s;k',s' ^^''('j-'k) ^s'' (^^k') r(Aa;) V'/p (Ak). 
^k,s;k',s' = Xi/fe (ek,s)* (ek',s')fc + (ek',s')J (ek,s)fc , (ek,ek') are unitary polarization vectors for signal and idler 
photons, ^^''(i^) — ~i^fiijj/2€{\L, s), if^ip (Ak) = cPr -tpip (r) exp (— jAk • r) and is the Laguerre-Gaussian pump 
field amplitude in E (r, z; t) = V/p (r) e*^'^^^-"^*)^!. r(Aw)=exp [iAlo {t - t/2)] sin (Aw t/2) / (Aw/2) is the time 
window function defining the Auj range given the interaction time r, Aoj = cjk + ujk — wp, Ak = k + k' — kp. Sym- 
metries associated with the crystalline medium influence Ak,s;k'.s' x 'tpip (^k) through x'^-*, x'^"* ^J^d the polarization 
vectors 0. 

It was derived in Ref. 2] that in order for a wave state of a one-photon field |V'(0)i to be an eigenfunction of Jz 
(or Lz for processes where Sz — 0, as will be assumed here) it has to obey Jz\tp{t))i = lTi\tjj{t))i. This directly gave 

I^W>i =E / d^kg{kp,kz,s;t)e''^a\k,s)\Q) . (3) 

In this condition one can see that the azimuthal phase should occur only as a phase term in the wave state and not 
in its amplitude. This guarantees that the magnitude square of the integrand presents complete rotational symmetry 
in (j) or, equivalently, presents a complete uncertainty in (f>, as required by Eq. ([1]). This condition has to be applied 
both to signal and idler photons. This implies that for perfect 0AM transfer from pump photons to SPDC photons, 
the probability amplitude i^s,s' (k,k') in Eq. ([2]) should contain the azimuthal angles for signal and idler only as 
phase terms e'""^ and e™"^ . As a linear superposition of wave functions is also a wave function, superpositions 
of expressions similar to Eq. ([2]) containing similar phase terms can also represent a valid solution. However, the 
superposition coefficients should not depend on the azimuthal angles. 

Starting from the accepted standard wave state for SPDC, one may expand the most relevant terms to compare 
with symmetry restrictions imposed by Eq. ([3]) and verify if a specific SPDC process may or may not transfer QAM 
to the down converted photons or even if this transfer can be partial. This answer should not be ambiguous and the 
result obtained should be general enough to allow direct comparison with specific experiments. 
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THE WAVE STATE AMPLITUDE F(ks,ki) 

The wave state amplitude i^(ks,ki) (see Eq. [2]) is the term that needs to be considered in detail. In principle, this 
amplitude contains a wealthy of information one could obtain from the wave state, from efficiency considerations to 
phase matching conditions. Of course, one has to assume complementary information about the crystalline medium 
including the light propagation conditions given by Fresnel equations. From the wave state amplitude the required 
conditions for phase matching emerge naturally. 

The term Ak,s:k',s' will be considered in a separate section and the time window term r(Aa;) admits easy inter- 
pretations. For example, for a CW laser where the interaction time r can be made large r(Aa;) — + tt6{Alu). The 
most involved term is the Fourier transform ipip{A\i). A careful consideration of this term allows one to derive main 
conclusions about the SPDC process in general. 



The Fourier transform i/'ip(Ak) 
The Fourier transform ?/';p(Ak) of ipip{A'k) can be written 



where 

ipip{r, (j), z) = 
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and = Wq [l + {z^/zj^}] , q{z) = z + izji, Wq — {2zR/kp) , zq gives the crystal center, = + j/^, and Ic is the 
crystal length along the propagation direction. Development of integral in r and (j) is straightforward although 
cumbersome: 



-0;p(Ak) = TV Alp 
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where pj. — Ak^ + Aky = + p'^ + 2pp'cos(0 — 0'), p — ksinO and p' — k' sinO'. 

The remaining z-integral can be solved under conditions favorable to experiments. Conditions zq = and Ic zji 
are most usual. This gives tan~^{z/{2zp)) and tan~^(zit:/z) 7r/2; therefore 



ijip{Ak)^-W+'2-%pL^]e'l^P+''^''"^^''-/^''-^]e-^/^e^^L'p (C)- 



izr- 
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, where i=^pl . (7) 
kp 



It is interesting to observe the presence of the phase tan^^(Afcj^/Afc2;) in ■0ip(Ak). It gives several possibilities for 
entanglement of signal and idler phases 4> and (j)' and is a signature of the complexity connecting signal and idler 
photons on the plane transverse to the propagation direction. The probability for unconstrained signal and idler 
occurrences |Fs^s'(k, k')p will be proportional to 



\^l;ip{A\^)\'^^'A-^\Aip\' ( ) e-^ e L^p (0 



ZR 



sm 



d^{A-2zRAk,+0 



(8) 



Eqs. (Ill) and ([8]) are fundamental for SPDC. Together with Ak,s;k',s' they determine probabilities of signal and idler 
occurrences and, by state superpositions with arbitrary phases, interferences. Phase matching are determined by the 
loci of maxima of these same equations. 
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PHASE MATCHING 

Phase matching conditions with variables A/c^ and ^ can be obtained from |^ip(Ak)p. In principle, polar and 
azimuthal angles that define maxima of |V'ip(Ak)p should be determined simultaneously but simplified solutions are 
frequently used according to the level of detail one requires. For example, since sinx/x in Eq. ([8|) is weakly dependent 
on ^, Eq. ([8|) can be treated as two independent parts for phase matching considerations. One part, in form of sinx/x 
gives the Akz range and the second or remaining part of the equation sets the width associated with the variable ^. 
Write |Vip(Ak)|2 ^ X where = (sinx/xf and = 7:H~'\Aip\^ {hzn/kpfe-^ ^ 4 {^f ■ /ion, can 

be considered non-negligible from Afc^^j,, = (4 + ^)/2zfl ~ to Ak^,,^^^ — {lc{4: + ^) — Attzr) /2l(,zii (for x — tt). For 
Ic <C 1, Akz — ±2tt/Ic- This sets the polar angles 9pm and d'p^^ for phase matching. Maxima for /t,ansv will describe 
azimuthal angles, (j) and </>', and 0pm and 0p„j, that are the polar angle values defined from Afc^^;,, to Afc^,,,^,;. These 
maxima are not hard to find for specific values of / and p. 

Just to exemplify a simplified use of Eq. ([5]) , one can calculate the probability for signal and idler occurrences in a 
specific Type II case of a uniaxial crystal where the pump beam carries a / = 4 0AM. The calculations are done for a 
BBO crystal with the crystal axis inclined by 6c with respect to the pump beam propagation direction. Two natural 
coordinate systems are involved, the crystal and the laboratory axes. While description of the linear and non-linear 
susceptibilities are usually given in the crystal reference system, the pump beam and the characteristic SPDC pattern 
desired define the laboratory system to be used. Appropriate coordinate transformations have to applied to provide 
correct answers. Starting with the longitudinal equation for phase matching, some additional information needed 
for its solution will be presented. These additional information will also be used for the transverse phase matching 
conditions. 



Longitudinal condition 

Expanding Akz in the longitudinal condition gives the interval where an appreciable contribution is found: Akz = 
±2'k/Ic. That is to say, for angles bounded by this condition, there is a good probability to find signal and idler 
photons. Let us re- write this condition in a more tight bound and expand Akz'. 

- kp + {kcos0 + k' cos0')~^ . (9) 

The wave vectors in the medium are k = (27r/A)n = {uj/c)n and k' ~ (27r/A')n' = (w'/c)n' where A and A' are vacuum 
wavelengths and the refractive indexes n and n' have to be found using Fresnel's equations for specific propagation 
directions and Sellmeier's equations for the principal refractive indexes. 



Fresnel's equations 



Fresnel's equations can provide the refractive index n of a uniaxial crystal along an arbitrary propagation direction 



specified by the wave vector k — (ujn/c) s where the unit vector is s = s^-x 
starting from Maxwell's Eqs. written in local variables: s-H = 0, D-s = 



Sy-y + SzZ. They are derived Q 
H X s = D/n , E X s = — H/rt, 
where V was written as V — > i{ujn/c)s. Eliminating H results Dj — [Ej — Sj (s • E)]. Considering an uniaxial 
crystal (e^j = euSij) and multiplying both members by sj one obtains £jjEj — n^Ej = —v? [sj (s • E)] that gives 

1= ' ' 



"'^s'j/{ejj — n^). Subtracting 1 (= s^) from both sides and writing = e, Fresnels's Eqs. arc obtained: 



— = 

1 1 " ' 



(10) 



where = sin cos 0, Sy — sinO sincf), Sz = cos 9 and e^^ 
principal axis and = Eg is the extraordinary dielectric constant 
Tz = e/ce and multiplying both members of Eq. (jlOp by (1 — rx)'^{l 



are the ordinary dielectric constants along the 
pn)) . Using the notation r^. — e/e^ 



Ty and 



one obtains 



(l-r^ril-rzY 



{l-rz){sl + sl) + {l~r^)sl 
(l-r,)(l-r,) 



. 



(11) 



FIG. 1: Decomposition of a wave vector k into transversal (p) and longitudinal (kz) components in the laboratory axis {x,y, z). 
A rotation around the j/-axis by 6^ separates the crystal principal axis from the pump beam propagation direction z. 

This equation defines that either {l — rx){l~rz) = 0, that gives either e = Co or e = eg (pure ordinary and extraordinary 
cases), or [(1 — r^) (s^ + s^) + (1 — r^) s^] = 0. This last equation can be written as 

1 si +31 si 

_ ^ y_ ^ ^ _ ,-^2) 

e Ce Co 

These are the possibilities for light propagation in uniaxial crystals. Either pure ordinary or extraordinary propagation 
or a mixture of ordinary and extraordinary light propagation. However, these equations define the refraction indexes 
in the crystal coordinate system. Obtaining the refractive indexes on the laboratory coordinate system demands that 
rotations are introduced corresponding to the geometrical situation adopted. 

Rotation matrices 

One should be aware that rotation angles do not obey a universal notation and sometimes references for crystal 
rotation angles may vary even from crystal to crystal in the literature. In particular, complementary angles may 
be referred to in a similar way. Fig. [T] shows a convenient laboratory coordinate system {x,y, z) obtained from the 
crystal axis by rotation of angle 9c around the y-axis. Usual crystal rotations to utilize better geometries to increase 
the efficiency of the down-conversion process consist of rotations about, say, one of the crystal secondary axis followed 
by a rotation on its principal axis. These matrix rotations can be written 

(cos <j!>c sin (pc \ / cos 6c — sin 6c \ 
-sin<j!>c cos<j!>e 0-01 . (13) 

1 / \ sin 6lc cos 6c / 

The rotation matrix connecting crystal and laboratory coordinate systems adopted in this work is defined by 
Eq. (jl3p . Of course, other commonly used rotation procedures exist. Rotation by two Euler angles, for example, are 
written not as Rzy but as Ryz, where the azimuthal rotation is applied first. These resulting rotations are not equal 
and they do not have a unique correspondence using just two Euler angles. Although what description to adopt is a 
matter of taste, it should be clearly stated to avoid misinterpretations about any result obtained. 

Refractive indexes 

A wave- vector unitary propagation vector k written in the laboratory coordinates (/>), k = 
(sin 6 cos (j), sin 6 sin cf), cos 9) is transformed to a vector in the crystal medium by s = Rzy ■ k. The resulting pump, 
signal and idler unitary wave vectors could then be plugged into Eq. (|12p to provide the corresponding refraction 
indexes in a uniaxial crystal such as the BBO. One obtains 



n = Tio 



(14) 



6 



np 



nf (cos 9c cos 6' + cos (j)' sin 9c sin 6'')^ 

np^e np^o 



-,/2 



(cos 6*' sin 9c — cos 0c cos (j)' sin 6*' ) + sin^ 9' sin^ 0' 



Up „ cos^ 



n D „ sin^ I 



(15) 



(16) 



It is easy to observe that propagation under eo polarizations do not present azimuthal symmetry. However, the overall 
symmetry for SPDC depend on other terms as well. 



Sellmeier equations 

The refractive indexes in the principal axes, np^e,np_OTno,nc,n'^, and n'g are usually obtained experimentally and 
represented by parametric equations representative of the microscopic physics involved. These equations, known as 
Sellmeier's equations, have a wide applicability due to its success to represent the refraction index of low absorption 
crystals as a real function of wavelength. They are based on the form of an oscillator response to an applied force 
or, as a microscopic theory for the response of bound electrons in a solid to an applied electric field (See [9], &2.3) 
r = eE/ [m (wg - ui^)] or P = Ner ~ x^^^^^ ~ {{n^ - l)/47r)i; and assumes the form =a + b/{\^ - c) - dA^. The 
constants a to d have to be experimentally determined for each crystal along ordinary and extraordinary propagation 
directions and at every desired temperature. An experimental fit to this phcnomcnological equation may give a 
very good numerical representation of the refractive indexes in a quite broad frequency range. The d term is a first 
corrective term of a possible series of even terms in A. The assumption of reality for the fields D and E and their 
causal connection imposes the even dependence on A for a real dielectric constant [lo| . This way, the first terms are 

nJXf = a - -J^ , and jUXf = e - —J- . (17) 

A^ — c A^ — Q 



Ref. [ll|, for example, present the experimental values a = 2.7405, b = 0.0184, c = 0.0179, d = 0.0155, e = 
2.3730, / = 0.0128, g = 0.0156, h = 0.0044 to represent BBO's refractive indexes along the crystal axes (A are given 
in /im in these equations) 



Longitudinal condition 



As the refractive indexes are now defined and the wave vectors k 
with adequate Sellmeier parameters, Eq. ([9]) 



kp 



{kcos9 + k' cos9') ~ — 

I'C 



{2n/X)n can be calculated for general angles 



27r 



can be solved. Due to the complexity dependence of the refractive indexes in the angles, this equation may not be 
solved analytically in a general case. Numerical solutions can be found and will determine the angles where the SPDC 
process is more efficient. For a Type II process in BBO, for example, Fig.[2]shows a plot obtained from the numerical 
solutions obtained. The obtained numerical solutions can be even parameterized for simplicity. This way, polar angles 
giving SPDC rings are closely represented by 



C(c 



7r) + v/c 



tt) +r] ] + i/exp [-lysin^ [(^s - 7r)/2]) cos [2{4>s - n)] 



and 



: arcsm 



C^cos^i + cos^ + rj^ +z/exp (— i^sin^ ((/ii/2)) cos2(/)i 



(18) 



(19) 



where C — 0.034, r] ~ 0.797, v ~ 0.0016, fi ~ 3.45. These equations give the thin black hues in Fig. [2l 

It has to be emphasized that these angles are described within the medium. Straightforward application of Snell's 
law gives the angles outside of the crystal. 



FIG. 2; Signal and idler rings obtained from numerical solutions and thin solid lines obtained from fitting the numerical 
equations with Eqs. ([TSj and ((19} . 




Transverse equations 



The equation /t„„„ = 7r^4~'|Aipp {lcZji/kp)'^e~^ (^)^ together with the obtained polar angles given by Eqs. 
(fT8| and (|19p define the loci of possibly entangled azimuthal angles that maximize Fs^^' (k, k'). Finding the analytical 
maxima for /transv is not a trivial task and will not be attempted here. However, for specific values of p and I this is 
usually a simple task. For example, Fig. [3] shows a plot of /transv/(7r^4~'| Ajpp (IcZu./kp)'^) for p = and I = A. It is 
easy to find the maximum at ^o4 — 15.9491. This defines the value pl{p = 0,1 = 4) = (fcp/zi?)^o4 that maximizes the 
signal and idler emissions. From the transverse equation 

pI = Akl + Akl = p2 + p'2 + 2pp' cos(^ - 0') - — Co4 , (20) 

and using the polar angle dependence given by Eqs. (HH]) and this equation can be solved for the azimuthal angles 
(f) and (f)' . Fig. [4] shows the numerical dependence found between these angles. While the condition 0' = + tt is an 
approximated one, some deviations exist that may be meaningful in some applications. These deviations are more 
apparent near azimuthal angles close to zero. The thin solid lines in Fig. [4] were obtained by fit, giving Eqs. pT|) . 

(t)\ = 3.1836 + O.3631e-°-46530^ + 0.9999(/) 

<^'_ = 3.0996 - O.3631e-°-46530' + 0.9999(/) . (21) 
The average azimuthal angle from these two equations is 

^' = ^^=0 + .. (22) 
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^'(rd) 




FIG. 4: Idler azimuthal angle versus signal azimuthal angle obtained numerically from Eq. (|20p . The thin solid lines were 
obtained from Eqs. 1211 




FIG. 5: Calculated transverse coincidence-count structures on s-ring with point detector on i-ring. All angles are lab angles but 
inside the crystal. The crystal is tilted with the laser beam at 9c = 49.7" ((f>c ~ 0) from the crystal c-axis. The laser wavelength 
is Xp — 351lA(ei = x)) and the principal refractive indexes are np,o = 1-707, np,e = 1.578, no = jIq = 1.665, = 1.548. 



Ref. pointed out that perfect transfer of 0AM from the pump to the SPDC photons demands azimuthal 
symmetry for |i^s,s' (k, k') ^ around the pump propagation direction (quantization axis). Lack of azimuthal symmetry 
causes partial transfer of 0AM. To illustrate this partial transfer. Fig. [5] shows transverse coincidence structures 
expected for degenerate non-coUinear Type-II SPDC in a BBO crystal (See for a Type-II experiment). Calculated 
structures represent the detection probability for signal and idler photons within small ALuALu'AdA6'A(j)A(j)' around 
phase matching conditions (Here sinx/x — > 1 and dk ~ {n/c)duj)): 



k,s:k' 



dwduj'dede 



(excluding existing geometric effects given by Ak,s;k' 



(23) 



The lack of azimuthal symmetry in Type-II SPDC is reflected on the coincidence structures (even with the neglect 
of |^k,s;k',s' P) that show highly asymmetric coincidence donut like pattern. Asymmetric structures should then be 
expected in Typc-II with 0AM. The Ak.s:k'.s' contribution (to be shown ahead) is non-negligible but presents no 
sharp variations for the structures. 
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indexes 


contractions 


mn 


11 22 33 23,32 31,13 12,21 


I 


1 2 3 4 5 6 



TABLE I: Contractions for tensor indexes. 



NON-LINEAR POLARIZABILITY 

The SPDC efficiency is directly proportional to the non-linear polarizability vector P = ^i,(k.s:k'.s')^*.cr with 

(2) 

components ^i,(k,s:k',s') given by the product of the tensor Xijk ^'^^ components of the unitary polarization vectors 

A,i)^,s-M',s') ^ xifk (ek,s)* (ek',s')fc + (ek',s')J (ek,s)fc ■ (24) 

The interaction energy Vj is given by the product of the laser field polarization and the non-linear polarizability vector, 
Vi = Ep.P. Up to know, all equations have been developed using laboratory coordinates. In order to keep the same 
reference system, all quantities in Ak,s;k',s' have to be referred to the laboratory coordinate system. Usually, the non- 
linear dielectric tensor is given in the crystal axes and a coordinate transformation to the laboratory axes is necessary. 
^k,s;k',s' involves the non-linear tensor x^^-* and the unitary polarization vectors e^ g. These constitutive elements 
will be considered in the next sections. The non-linear tensor can be written in the laboratory coordinate system and 

(2) 

multiplied by the components of the unitary polarization vector in the same reference system or, alternately, Xijk ^'^^ 
the unitary polarization vectors can be written in the crystal reference system and the resulting vector component 
rotated to the laboratory coordinate system. This last method will be followed. 

(2) . . . . (2) 

The susceptibility tensor Xqmn for the class of uniaxial crystals can be written in a contracted form as Xqmn ~^ 
xfi = 2dq/ as indicated in Table HI This way. 



^ dii di2 di3 di4 di5 die 

^21 ^22 ^23 d25 \ (25) 

\ ^31 ^32 ^33 1^34 £^35 '^36 



To calculate the unitary polarization vectors, one may write the electric displacement D in the crystal principal 
axis, where = CuSij results Di = eaSuEj — [E^ — SiSkEk] or 

[nl - n'^ (l - si)] Ex + n^SxSyEy + n'^Sx-s^.Ez = 
n'^Sx-SyEx + [uy - (l - s^)] Ey + n^SyS^E^ = 

n^s^s.Ex + n'^SyS.Ey + [nj - n^l ~ s^)] E, ^ . (26) 

This set of Eqs. give the possible electric field amplitudes for specific propagation directions and refractive indexes n. 
Without developing general solutions, one may look for a solution for the signal photons propagating with ordinary 
refraction index tIq. Writing = riy = rio and = rig for the idler's refractive indexes, the electric field components 
in the medium are obtained and normalized resulting in 

Go, cr =(— singer, cos 6'cr, 0) . (27) 

To calculate the unitary polarization vector for extraordinary propagation, one may observe that in the crystal 
medium the Poynting vector form a set of orthogonal axes with the electric field and the magnetic field. At the same 
time, the electric displacement may not be along with the electric field but is normal to the propagation vector in the 
medium. Eq. (j26p can be particularized for this case as well. However, for simplicity, one may look for field solutions 
that give a unitary polarization vector for the extraordinary propagation e^^ orthogonal both to and to its unitary 
propagation vector s' = s^x -I- sj^y -|- s'^z. This gives 

^'ecr = (~ COS 9'^^. COS (j)'^^. , — COS 6'^^. SlU (j)'^^ , siu cos ) . (28) 
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Non-linear polarizabilities in the crystal reference system 



Calculation of the non-linear polarizability components ^i,(k,s;k',s') demands specification of 
(ek,s)* (ek',s')fc + (ek',s')! (ek,s)s 



for all SPDC possible cases. The question mark in Eq. (p9| represents the 
need for these choices. The possible cases (0,0), (e,e) , (o, e) and (o, e) are detailed in the Appendix. 



(k,s;k',s') 



A 



2,(k,, 



\ ^3,(k,s;k',s') 



^ dll di2 
^21 C?22 
\ rfsi ^32 




/ eie'i \ 

6262 
6363 

6263 + 6362 
6163 + e^e'^ 
V 6162 -I- 6261 / 



(29) 



Using the appropriate polarization cases the signal and idler non-linear polarizability Poo, Pee, Peo and Poe are 
calculated (see Appendix). These non- linear polarizabilities cover all possibilities of light excitation in uniaxial crystals. 



Non-linear polarizabilities in the laboratory reference system 



For SPDC it is usual to have the crystal rotated by convenient angles referred to the laboratory reference system 
X, y, X, 9, 4> instead of Xcr, 2/cr, Xcr, dcr, 4>cr- A laser with amplitude Ep polarized along direction e will define the light 
matter interaction to be analyzed. In this work it is chosen to have the non-linear polarizability P rotated to the 
laboratory reference system; this transforms P Pjab- 

Eqs. ()43|44|46|48|) can be written in the laboratory coordinate system under rotation given by R~y] PaPiiab — 
R^y.Vap, {aP — oo,ee,eo,oe). For an example, one can choose a laser amplitude Ep = Ep'i and obtain for o, e 
polarizability 

Vi^oe = —EpSi ■ (i?^J^.Poe) = cos 6*0 COS (/>c (—^14 COS ^c- sui 9'^j. + di2 COS 9'^^. cos (f)cr sin 0'^^ 

—dll cos9'^^ cos0or singer + di5 sin 6*^,. sini/icr + c?i6 (cos 6*^,. cosc/)^,. cosher — cosO'^j. smcj)'^^ singer)) 
— cos 0c sin 0c (—^24 cosher sin 6*^^ + ^22 cos^^r cosher sin0'^^ — ^21 cos6'cr cosher singer 
+d25 singer singer + (^26 (cosher COS 0'^^ COS (per ~ COS 9'^^. sin 0'^^ singer)) 
+ sin 9c (~ (<i34 COS ipcr sin6'^^) + ^32 cos 9'^^ cos (per sin^'^^ — dsi 0036*^^ cos ip'^^ sin (per 
+d35 sin6'^^ singer + ^36 (cos0^^ cos^'^^ cos(pcr — 0036*^^ sin(p'^^ silKpcr)) 

For a crystal where the dominant coefficients are dn, ^22 and dis (e.g. BBO), V/^oe simplifies to 

Vi^oe = — {d22 COS 9c COS 9'^^ COS (per sin (pc sin (j)'^^) + cos 9c cos (pc {—dll cos 9'^^ cos (p'cr sin (per + rfi5 sin 9'cj. sin (per) ■ (31) 

Replacing 9'^^ ^ 9er (per ^ <^cr in Eq. ([3T|) gives Vj^eo for BBO. 



Crystal to laboratory reference system 



If one wishes to have Vf^oe (or any other interaction energy term) written under laboratory angles {9,(p), the con- 
nection between angles in these two systems have to be found. Given a unitary vector v = (sin 9 cos (p, sin 9 sin 0, cos 9) 
in the laboratory reference system, the rotation R-^y brings it to the crystal coordinate system (or the inverse rotation 
to move the reference system). Assuming that vector components in the crystal reference system can be written as 
V — Wer — (singer cos 0cr, singer sin 0cr, cos 6*^), a connection between the medium and the laboratory angles can be 
found. This way 



9er — arccos [cos 9e cos 9 + cos (p sin 9e sin 9\ 



\/n\ ^ ni 



■\/(— 1 + cos 9e cos 9 \- cos (p sin 9e sin 9) (1 -|- cos 9c cos 9 -|- cos (p sin 9c sin ( 



(32) 
(33) 
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where 

nl = — 1 + cos 0^ sin sin 6*^ + cos (cosS^ + coS(?i>^ sin6'^ sinc^c^) (34) 

n2 ~ {cos 9 sin 9c sin 4>c + cos 0c sin0 sin0)^ (35) 

n3 = 2cos6'c cos (pc cos0 sin0 (cos 6' cos 0c sin6'c — sin6' sin 0c sin0) . (36) 

Replacing in Vj^oe the appropriate angle given by the 9cr or 0cr above, Vi^oe is determined in terms of the laboratory 
angles. The procedure is straightforward but the resulting expression is quite long. 



A method for determination of non-linear coefficients 



Knowing the precise rotation angles between the crystal reference system, the laboratory system and phase matching 
angles, it was shown that all elements to determine the signal and idler scattering probabilities (wave state amplitude, 
Eq. ([7])) can be obtained in detail. From these information, transverse coincidence-structures can be calculated (see 
0) and by integration over the angle variables for signal (or idler) the angular dependence of the idler (or signal) 
intensity can be obtained. Aside from obtaining a detailed picture of SPDC process, it should be pointed out that use 
of these equations allows one to obtain directly through simple intensity measurements [l3| the non-linear coefficients 
dij. This can be done by fit of the existing experimental angle dependence to the theoretical elements presented, 
including numerical integrations indicated. Although this analysis is not the object of this work, it is interesting 
to observe that considering just an integration over the nonlinear polarizability, neglecting the contribution from 
'0ip(Ak), leads to azimuthal asymmetries. For example, on the transverse plane {x,y), normal to the pump laser, 

Poe,x,y = [cos 0c (^15 COS 0c siu 6*^^ siu 0cr — COS ((^22 COS 0cr Siu 0c siu 0'^^ + du COS 0c COS 0^^ sin0cr)) , 

d22 cos6'cr COS0C COS 0cr siu 0^^ + (— {du COS 6'cr cos0cr) +^15 sinfl^^) sin0c sin0cr,O] , (37) 



27r \ 2 

Poe,x,yd(j)'cr I = 4d^57r^ siu^ siu^ 0cr (cos^ 9c cos^ (j)c + sin^ 0c) . (38) 



The (simplified) azimuthal dependence indicated by Eq. (|38p shows a variable intensity for the SPDC rings (See 
experimental result in [l3|). 



Modulation by external fields 



Extending our knowledge about the microscopic behavior of the light-matter interactions in SPDC would help 
us to examine other possibilities to use entangled photon pairs. External generalized fields can be added to allow 
modulations to be applied to these systems. Pressure, electric and magnetic fields are natural candidates to exert 
different modifications on the optically nonlinear medium. The virtual character of the SPDC process do not allow 
direct access to energy levels involved (At — > leads to Ai? oo) but, nevertheless, it does not exclude one to 
observe important effects related to these virtual processes. For example, a magnetic field may modify electronic 
levels with detectable effects on SPDC (Type I or Type II). Induced defects can also be used to probe for local 
symmetry variations in crystals 0. Several tools can be used to explore this fundamental problem of 0AM transfer 
by a non-linear medium and may lead to a better understanding of the underlying microscopic physics. Future 
quantum applications of 0AM entanglement, such as quantum computation or teleportation, may depend on a deep 
understanding of these 0AM transfer process to achieve a very efficient use and control of quantum entanglements. 



CONCLUSIONS 



Explicit calculation of the equations determining SPDC processes when 0AM is involved were provided for crystals 
of uniaxial symmetry. The light-matter non-linear polarizability components ^i,(k,s;k',s') were calculated giving the 
complete angular dependence in the crystal reference system as well as in the laboratory system. The Fourier 
transform '0ip(Ak) was analytically calculated and approximations used were discussed. Phase matching conditions 
were obtained and it was shown that spatial transverse coincidence structures can be calculated. In Ref. ^] it was 



12 



shown that whenever |F(ks,ki)p (see Eq. ([2])) lacks azimuthal symmetry, an expansion of F(ks,ki) in terms of the 
azimuthal angles reveals that -despite the energy conservation condition- the initial orbital angular momentum I 
connected with a photon in the incoming mode may not be completely transmitted to the SPDC pair. The transfer 
may be partial, with the signal and idler carrying the 0AM value I' ^ I. The obtained equations for ^i,(k,s;k',s') and 
'0ip(Ak) allows one to make these expansions to study specific cases. It is expected that the explicit treatment of the 
probability amplitude given in this work may allow further developments in the study of quantum images in SPDC 
process involving 0AM. Analytical tools are then provided to indicate whether a specific SPDC process may or may 
not transfer 0AM to the conjugate photons. Straightforward extensions of this work can also be done such as to 
obtain output profiles in second harmonic generation where one or both of the input beams are in 0AM states. The 
symmetry properties of the light-matter interaction will be revealed by the up-converted beam. 
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APPENDIX 



Expansion of the vector products give the cases (o, o), (e, e) , (o, e) and (o, e): 



\ 



eiei 

6262 

6363 
6263 -I- 6362 
6163 -I- e36'i 

y 6162 -I- 626i j 



60,160,1 
6o,26o_2 
60,360,3 
60,260^3 -I- 60,360,2 
60,160,3 -I- 60,360,1 
\ 60,160,2 + 60,260,1 / 



y — cos ( 



Sm (per sin (por 
COS ^cr COS ^cr 





sin (;/>or ~ sin (pc.r cos ^^r / 



(39) 



/ 



6161 
6262 
6363 
6263 -I- 6362 
6163 -I- e3e'i 

y 6162 -I- 626i ) 



I 6e,l6e,i 
6e,26e,2 
6e,36e,3 
6e,26^,3 -I- 6e,36e,2 
6e,l6e,3 -|- 6^,36^,1 
V 6e,l6e,2 -|- 6e,26e,i ) 



^ COS dcr COS Q'^^ COS ^cr COS ^^r 

cos Oct cos Q'^^ sin t^cr Sin (;/>or 

sin QcT sin Q'^^ 
— sin Oct cos Q'^^ sin t^'^^ — cos Qcr sin t^cr sin < 



— cos 6 
y cos QcT 



cos 



cos ^cr sin ^ sin Qcr cos COS 0^, 

cos (/icr sin (/i^r + sin ^cr cos ^or ) / 



(40) 



6161 
6262 
6363 
6263 -I- 6362 
6163 -I- 6361 

y 616^, 626'i ) ^ ^ 



I 6e, 160,1 \ 

6e,26o,2 
6e,36o,3 
6e,26o,3 -|- 6e,36o,2 
6e, 160,3 -I- 6e,36o,l 
V 6e,l6o,2 -I- 6e,26o,l / 



— COS Qcr cos <^cr Sin (^'cr ^ 
cos Qcr sin 0cr COS (\)i.r 


— sin Qcr cos i^'cT 
sin Qcr sin 

y cos Qcr COS 0cr COS 0^.^ — COS Qcr Sin 0cr Sin j 



(41) 



6161 
6262 
6363 

6263 4- 6362 
6163 -f 6361 

y 6162 -I- 6261 j 



( 60, l6e,i \ 

6o,26'e,2 
6o,36e,3 
6o,26'e,3 -|- 60,365,2 
6o,l6e,3 -I- 60,365,1 
y 6o,l6'e,2 + 6o,26'e,i j 



\ 



— sin (\>cr COS Q'cr COS (\>cr ^ 
COS 0cr COS Q'cr Sin ^^r 


— COS t^cr sin Q'cr 
sin (jicr sin Q'cr 

COS 0cr COS Q'cr COS l\>cr " sin (l>cr COS Q'cr sin j 



For the coUinear and degenerate propagation, these vectors give the particular cases described in Rcf. [12 
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The nonlinear polarizabilities are obtained in a straightforward way, giving 

Poo = [sin0'c^ (dii singer - die cosher) + cos^or (^^12 cosher - die singer) , 

Sint/)^^ (^21 sin (jjcr - ^26 cosher) + COScj)'^^ (^22 COS (jjcr - ^26 Sin (j)cr) , 

sin0^^ (dai singer - ^36 cosher) + coscf)'^^ {d^2 cosher - ^36 Sin(t>cr)] , (43) 

-f*ee \-Pee.xcr 1 Pee^ycri -^ee.zcr^ i where (44) 

Pee,xcr — sinO'^^. ((^13 sin — COS 6'cr (c'ls cos4>cr + c'i4 sin^or)) + COS 6*^^ (— sin 6'cr (^15 COS + (^14 sin^^^) 

+ cosher {dll COS(/)or cosher + ^16 COS (/)cr Sin 0^,. + dig COS 0^,. sin 0cr + di2 Sin0^^ sin (f>cr)) , 
Pee,ycr = sin 6*^,. (o?23 sin 6'cr " COS 6'cr (^25 COS 0cr + 6^24 sin 0cr)) + COS S^r (~ sin 0cr (^25 COS 0cr + 1^24 Sin^^r) + 

cosher (^21 COS(j)'^^ COS 0cr + ^26 COS (j)cr sin + 1^26 COS sin 0cr + C?22 sin sin (j)cr)) , 
Pee,zcr = sin 6'cr (^33 sin 6'cr " COS 6'cr (^35 COS 0cr + ^34 sin^cr)) + COS 6'cr (" sin 6'cr (^35 COS ^c^ + ^34 sin (f)'^,,) + 

COS 9cr {dsi cosher COS (/)cr + ^36 COS (/)cr sin 0cr + die COS (t)'^^ sin (/)cr + ^32 sin(t)'^^ singer)) • (45) 



Peo \Peo.xcr 1 Peo.ycri Peo^zcr\ ; where (46) 

Peo,xcr = sin 6'cr (-di4 COS (/t'c^ + cii5 sin (/)cr) 

+ COS 6'cr (cos (/)cr (^16 COS (j)cr + ^12 sin (j)cr) — sin (j)'^^ [dn COS (j)cr + die sin 4>cr)) , 
Peo,ycr = sin 6'cr (-0^24 COS (f)'^,, + d25 sin (f)'^,,) 

+ COS 6'cr (cosher {^26 COS (j)cr + ^22 sin (j)cr) - sin 0^,. (^21 COS (?!>cr + ^26 sin (j)cr)) , 
Peo,zcT = sin 6'cr ("^34 COS (/)cr + '^35 sin 0cr) 

+ COS 6'cr (cosher (1^36 COS (t>cr + ^32 sin 0cr) " Sin 0cr (c^31 COS 0cr + ^36 sin (/)cr)) • (47) 

Poe — \Poe^xcr 1 Poe,ycr ^ Poe,zcr\ : where (48) 

Poe,xcr ^ sin 6'cr (-^14 COS (/)cr + ^15 sin 4>cr) 

+ COS 6'cr (cos 0cr {die COS(/)cr + di2 singer) ~ ('^ll COS 0cr + ^16 sincj)'^^) sin (/)cr) , 
Poe,ycr = sin 6'cr (-^24 COS (/)cr + ^25 sin 0cr) 

+ COS 6'cr (cos (^cr (^26 COS (j)'^^ + 1^22 Sin0^r) ~ (^21 COS 0cr + ^26 sin (/)cr) sin (^cr) , 



Poe,zcr = sin 6'cr (-(^34 COS (/)cr + ^35 sin (/)cr) 

+ COS 6'cr (cosher (^36 COS (j)'^^ + ^32 sin(/)cr) — (^31 COS (/)cr + ^36 sin ( 



cr 



.) sin(/)cr) . (49) 
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